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^ Abstract 

Let R be a real closed field (e.g. the field of real numbers) and ^ C R" be a semi-algebraic 

Yl set defined as the set of points in R" satisfying a system of s equalities and inequalities of 

I— ^ multivariate polynomials in n variables, of degree at most D, with coefficients in an ordered 
ring Z contained in R. 

►^ We consider the problem of computing the real dimension, d, oi ,y. The real dimension is 

QQ the first topological invariant of interest; it measures the number of degrees of freedom available 

OQ to move in the set. Thus, computing the real dimension is one of the most important and 

^\ fundamental problems in computational real algebraic geometry. 

' I The problem is NPR-complete in the Blum-Shub-Smale model of computation. The current 

■^^ algorithms (probabilistic or deterministic) for computing the real dimension have complexity 

O (si:))0(d(n-d))^ ^^^^ becomes (s £))'='("') in the worst-case. 

, The existence of a probabilistic or deterministic algorithm for computing the real dimension 

• • with single exponential complexity with a factor better than O(n^) in the exponent in the 

, J^ worst-case, is a longstanding open problem. 

S^ We provide a positive answer to this problem by introducing a probabilistic algorithm for 

?H computing the real dimension of a semi-algebraic set with complexity {s D)'^^"K 



1 Introduction 

A semi-algebraic set S^ C R", where R is a real closed field, is defined as the set of points in R" 
satisfying a Boolean formula whose atoms are polynomial equalities and inequalities. 

Computational real algebraic or semi- algebraic geometry is the study of effective algorithms for 
computing with semi- algebraic sets. Besides being a fascinating and important research area on its 
own, it is also one of the cornerstones of theoretical computer science. 

Many important results rely on the foundations of real algebraic geometry. Let us mention non- 
linear computational geometry [20, 44, 45], the recent breakthroughs in combinatorial geometry 
on the discrete version of Kakeya problem [30, 34, 47], and the new algorithms for non- negative 
matrix factorization [1, 38] based on testing the emptiness of semi-algebraic sets. Last but not 
least, we emphasize the intrinsic connection between computational real algebraic geometry and 
game theory, especially stochastic games [17, 21, 31, 39, 48]. 

Typical computational challenges in real algebraic geometry are algorithms for deciding the 
emptiness and/or computing at least one point at each semi-algebraically connected component of 
a semi-algebraic set [11, 13, 26, 29, 41, 42], algorithms to perform geometric operations such as 
projection (this operation is tightly coupled with quantifier elimination) [11, 13, 28, 41], answering 
connectivity queries (roadmaps) [10, 20, 32], computing the real dimension of a semi- algebraic set 
[35, 36, 50] or computing more sophisticated topological information, such as the number of semi- 
algebraically connected components, the Euler-Poincare characteristic, Betti numbers [8, 9, 12, 13]. 

Denote by s the number of polynomials involved in the description of a semi- algebraic set, by 
n the number of variables, and by D the maximum of the degrees of these polynomials. We can 
solve almost all the problems in computational real algebraic geometry using the generic approach 
of cylindrical algebraic decomposition [24] albeit in double exponential time, (s D)^ . Even 
though huge effort has been invested the last 25 years to derive algorithms with single exponential 
complexity w.r.t. the number of variables, there are problems that are still missing an algorithm with 
such a complexity bound. Moreover, even in the case where the complexity is single exponential, 
the exponent is not always 0{n). 

Let us emphasize that improving the exponents in the complexity bounds of algorithms in com- 
putational real algebraic geometry is not only a theoretical challenge. It introduces new algebraic 
and geometric techniques that find applications in more general domains, and eventually leads 
to efficient implementations for real-world problems. For instance, the first improvement of the 
long-standing 0{n'^) exponent in the complexity bound of Canny's probabilistic algorithm [20] to 
0{n'^''^) is based on a new geometric connectivity result that introduced the use of a baby steps 
giant steps algorithmic technique in this problem [43]. 

On the other hand, the problem of computing the real dimension lacks, up to now, an algorithm 
with single exponential complexity and exponent 0{n). 

Problem statement and state-of-the-art. In this paper we address the problem of computing 
the real dimension of a semi- algebraic set. The following definition is in order: 

Definition 1. [13, Section 5.3] Let S^ he a semi-algehraic-set of R" , where R is a real closed field. 
The real dimension of S is the largest integer d such that there exists an injective semi-algebraic 
map from (0, 1) to S^ . By definition the dimension of the empty set is —1. 

The best known complexity bound, in the worst case, for computing the real dimension of a semi- 
algebraic set is due to Koiran and it is (s D)^"' ' [36], where s is the number of polynomials used to 



describe the semi-algebraic set. It is based on quantifier elimination techniques, see [13, Alg. 14.10] 
and references therein. A partial improvement of the O(n^) in the exponent is due to Vorobjov 
[50]. He presented an algorithm with complexity (sL))*^('^("~'^)), where d is the real dimension. 
This bound is output sensitive, and when d is a constant, then it becomes (sL>) '"^ Basu, Pollack, 
and Roy [14] slightly improved the result of Vorobjov, based on [22]. They presented a complexity 
bound that depends on whether d > n/2 oi d < n/2. This result has a better dependence on the 
number of polynomials, s, than the one of Vorobjov [50]. 

On the other hand, it is well understood that we can compute the (Krull) dimension of an 
algebraic variety over algebraically closed fields in time D^"'' [23, 26], see also [35]. In the alge- 
braically closed field case, we can consider a sufficiently generic, random, collection of d hyper planes, 
for 1 < d < n, and check whether their intersection with the algebraic set under consideration is 
finite. The largest d where this is achieved imposes that the Krull dimension of the algebraic set 
is d. However, when we are interested in computing the real dimension of a real algebraic, or 
semi- algebraic set, this strategy is not applicable. 

It is of great interest to know if the problem of computing the dimension admits the same 
complexity bound in the real case and in the algebraically closed case. Quoting Koiran [36] "The 
main open problem is whether DIMjr [the real dimension] can be solved in time {s Z))*^(") ". Vorobjov 
[50] also mentions that "For a real variety V existence of a probabilistic dimension algorithm, with 
complexity bound {s D)^^^' is an open problem". 

The purpose of the present work is to provide a positive answer to this open problem. 

Besides the intrinsic mathematical interest for an improved algorithm for computing the real 
dimension, improvement of the complexity bound has important consequences. Some algorithms 
in computational real algebraic geometry consider the real dimension of a semi- algebraic set as 
part of their input, e.g. [13, Theorem. 13.37]. Let us also mention the recent bounds in [6, 7] on 
the number of semi-algebraically connected components of a semi- algebraic set that depend on the 
real dimension of some real algebraic set. Moreover, effective algorithms for computing the real 
dimension are needed to estimate efficiently the degrees of freedom in robotic mechanisms (see e.g. 
[33, 40]). 

The problem is also very important from the complexity theory point of view, as Koiran [36] 
proved that it is NPjj-complete in the Blum-Shub-Smale computation model [16] . 

Our results. We present an efficient algorithm for computing the real dimension of a semi- 
algebraic set. 

Our algorithm reduces the unbounded case to the bounded one, using a standard technique of 
computational real algebraic geometry introduced in [11]. 

Previous approaches for computing the real dimension of a semi- algebraic set 5^ rely on finding 
the largest integer d for which there exists a d-dimensional linear subspace, such that the projection 
of S^ on this subspace has dimension d. To do so they use quantifier elimination. Therefore, the 
exponent in the complexity is the dimension, d, multiplied by the number of quantified variables, 
n — d. 

The algorithm that we present, instead of projecting the semi-algebraic set under consideration, 
exploits geometric properties of fundamental objects of algebraic geometry, that is polar varieties. 
Roughly speaking, polar varieties are the critical loci of projections, e.g. [3, 4, 5] and references 
therein. More precisely, we are able to prove the following: Let V be algebraic set defined by the 
polynomial equalities of the input and U the open semi-algebraic set defined by the polynomial 



inequalities of the input. Then, up to a generic change of coordinates, d is the largest integer such 
that J/" is equal to the intersection of U and the limit of the critical locus of the d + 1 polar variety 
V, after we perturbed it symbolically. This way the computation of the real dimension reduces to 
finding the largest integer d with this property. 

The algorithm is probabilistic since we perform a random linear change of coordinates in the 
beginning. If we work over the integers, then we denote by r the maximum bit size of the coefficients 
after the linear change of coordinates. 

Our main result is encapsulated in the following theorem: 

Theorem 1. Let F = {fu ■ ■ ■ , fp) € Z[Xi, . . . ,Xn]P, G = (51, . . . , <7s) € Z[Xi, . . . ,X„]^ where the 
degree of each fi, resp. gj, is at most D. There exists a probabilistic algorithm for computing the 
real dim,ension of the semi-algebraic set defined by 

/i = • • • = /p = 0, 51 > 0, . . . , 5s > , 

in (s D)^^""' operations in Z. 

// Z = Z, then the Boolean complexity of the algorithm r (s D)^^' . 

To the best of our knowledge this is the first algorithm for computing the real dimension of a 
real algebraic or semi- algebraic set within this complexity bound. 

Organization of the paper. The rest of the paper is structured as follows: In the next Section 
we present the necessary preliminaries from real algebraic geometry. In Section 3 we present the 
genericity properties that the semi-algebraic set under consideration should satisfy. We prove that 
a semi-algebraic set can satisfy these properties if we apply a random linear change of coordinates. 
Section 4 presents the geometric result that is the crux of the matter of our algorithm. Finally, 
in Section 5 we present the algorithm ComputeRealDimension for computing the real dimension, its 
various subroutines, the proof of correctness and the complexity analysis. 

2 Preliminaries 

In this section, we introduce some basic notions and some notations that are used throughout the 
paper. 

As sketched in the introduction we will introduce some infinitesimals to deform real algebraic 
sets. This will lead us to consider various ground fields and semi- algebraic sets defined over these 
fields. We refer the reader to [13, Chapter 2] for a more detailed exposition of these notions on 
real fields, real closed fields, infinitesimals and semi-algebraic sets. We will also use basic notions 
coming from algebraic geometry since we will use the knowledge of the dimension of some algebraic 
sets to deduce the real dimension of semi- algebraic sets under study. For a more detailed exposition 
of these notions, we refer the reader to [46, Chapter 1]. The section finishes with some notions on 
critical points and polar varieties that are extensively used in the sequel. 

Ground fields. Let Q be a real field, R be a real closed field and C be the algebraic closure of 
R. We consider a field K containing Q (e.g. R or C) and let e be an infinitesimal. In the sequel, 
]K(e) stands for the Puiseux series field. We say that z = ^j>j ojE*''^ G ]K(e) is bounded over K if 
and only if zq > 0. We say that z = (zi, . . . , Zn) £ IK(e)" is bounded over K if each Zi is bounded 



over K. Given a bounded element z S ]K(e) , we denote by lini£_i.o z the number ao G IK. Given a 
bounded element z G IK(e)", we denote by lime_>o z the point (lim£_j.o(2;i), . . . ,lime_5.o(2;n)) G ]K". 
Given a subset ^ C IK(e)'", we denote by lim£_j.o(^) the set {lime_^o(^) \ z £ A and z is bounded}. 
Given a semi- algebraic (resp. constructible) set ^ C R" (resp. A C C") defined by a quantifier- free 
formula $ with polynomials in R[Xi, . . . ,X„], we denote by ext(^, R(e)) (resp. ext{A, C{e))) the 
set of solutions of $ in R(e)" (resp. C(e)"). 

In the sequel, we will work with n-variate polynomials with coefficients in Q, Q[C] and Q[e,C] 
or e and C are infinitesimals with < e < (^. Sign conditions on finite families of polynomials with 
coefficients in Q define semi- algebraic sets in R", those with coefficients in Q[^] (resp. Q[e, Q) define 
semi- algebraic sets in R(C)" (resp. R(0(e)"'). 

Basic definitions on algebraic sets. Let IK stand for an algebraic closure of K. We consider 
algebraic sets in K" defined by polynomial equations with coefficients in IC. A Zariski open set is 
a set whose complementary is an algebraic set. A constructible set in K" is the set of common 
solutions of a system of a finite number of n-variate polynomial equations and inequalities with 
coefficients in K. 

Let V C K" be an algebraic set defined by polynomial equations in ]K[Xi, . . . , Xn]- 

We will consider the dimension of V, referring to its Krull dimension (see e.g. [25]). This notion 
of dimension coincides with other notions coming from differential or algebraic geometry (see e.g. 
[25, Part II]). Roughly speaking, it is the number of generic hyperplanes such that their intersection 
with y is a finite set of points. The Krull dimension of a constructible set is the Krull dimension 
of its Zariski closure. If W is another algebraic set and V C W , then the Krull dimension of V is 
less than or equal to the Krull dimension of W . 

The algebraic set V is said to be irreducible if it cannot be decomposed as the union of two 
algebraic sets defined by polynomial equations with coefficients in K. If V is not irreducible it 
can be uniquely decomposed as a finite union of irreducible algebraic sets; these sets are called the 
irreducible components of V . 

When all the irreducible components of V have the same Krull dimension, we say that V is 
equidimensional. The ideal associated to V is the set of polynomials with coefficients in IK which 
vanish on V . There exists a finite family of polynomials which generate it in IK[Xi, . . . , X„]; let us 
denote it by /i,...,/p. 

Assume that V is equidimensional of Krull dimension d. A point x G 1/ is called regular (or 

smooth) if the Jacobian matrix 

r dh dfi 

dfp dfp 

has rank n — d at x. The kernel of the above Jacobian matrix at x is the tangent space to V at 
x; we denote it by TxT^. Points in V that are not regular are said to be singular. An algebraic set 
with no singular points is smooth. 

Semi-algebraic sets and algebraic sets. Let C be the algebraic closure of R, and S* C R*^ be 

a semi- algebraic set. The smallest algebraic set containing S is called the Zariski closure of S. It 
is well-known that the Krull dimension of the Zariski closure of S equals the real dimension of S 



(see e.g. [19, Proposition 2.8.2]). In particular if W is an algebraic set that contains S, one can 
conclude that the real dimension of S is less than or equal to the Krull dimension of W. 

Let 5 C R"' and S" C R". Consider a semi- algebraic map (/? : 5 — >• S" and R' be a real closed field 
containing R. We will consider the extension of (p to R', denoted by ext(c^, R'), as the semi- algebraic 
function ext(5, R') — )• ext(5", R') whose graph is the extension of the graph of ip to 5'. 

If 5 C 5", then by Definition 1 the real dimension of S is less than or equal to the real dimension 
ofS'. 

If X G 5, we say that x is a smooth point of S if it is smooth in the Zariski closure of S; a 
semi- algebraic set is said to be smooth if it is contained in the set of regular points of its Zariski 
closure (this is direct a consequence of [19, Definition 3.3.4]). 

We also consider canonical projections ttj : {xi, . . . , Xn) i— ^ (2:1, . . . , Xj). 

Change of variables. Let K be a field containing Q and K be its algebraic closure. Consider 
/ G ]K[Xi, . . . ,X„] and V C K" be the set of roots of / in K and A G GL„(Q). We denote by 
/ the polynomial /(AX) and by V C K". In other words, V is the image of V by the map 
x — )• A~^x. 

Similarly, we will also consider change of variables on semi- algebraic sets. If S" is a semi- algebraic 
set in R" and A G GL„(Q), then S^ denotes the image of S by the map x — t- A^^x. 

Critical points and polar varieties. Let fi, . . . , fp be polynomials in ]K[Xi, . . . , X„] and V C 
C" be the algebraic set defined by /i = • • • = /p = 0. For 1 < i < n — p, we define the polar variety 
crit(7rj, V) as the set of points in V at which all p-minors of the truncated Jacobian matrix 

8/1 ... dh_ 

dXi^i dXn 

dfp . . . dfp 

dXi^i dXn 

vanish. By convention for i > n — p + 1, we set crit(7rj, V) = V and for i = we set crit(7r,, V) = 0. 

We refer to [5] for a detailed study of polar varieties. Below, we recall some basic results that 
are used in the paper. 

Given a polynomial / G ]IC[Xi, . . . , X„\ and V C K" defined by / = 0, we denote by crit(7rj, V) C 
K" defined by 

; = ^^=... = ^=„. 

dXi+i dXn 

Assuming / is square-free, the set of singular points of V is defined by the vanishing of F and all its 
partial derivatives. The set crit(7rj,y) is called polar variety associated to ttj. When V is smooth, 
crit(7rj,y) is the set of critical points of the restriction of ttj to V (i.e. the set of regular points 
X G y such that 'Ki{Tyy) has dimension < i — 1). When V is not smooth crit(7rj, V) is the union 
of the singular points of V and the critical points of the restriction of ttj to V . 

Let us also mention that for i = 1, crit(7ri, V) contains the local minimizers and maximizers of 
the restriction of vri to F n R". 



3 Genericity properties 

In this paper, a property is called generic (in some suitable parameter space) if it holds in a 
non-empty Zariski open subset of the parameter space under consideration. 

Definition 2. Let f G Q((")[Xi, . . . , X„], V C C(C)"' be the algebraic variety defined by f = 0, and 
Ve C C(C)(e)"' be the algebraic variety defined by f — e = 0. We say that f satisfies property N if 
the following conditions hold: 

Ni for all 1 < i < n, crit(7rj, V^) is either empty or is smooth and equidimensional of Krull dimen- 
sion i — 1; 

N2 for a// X E y n R(C)"; tt^ (^d(x)) H (V Cl R(0") is finite where d is larger than or equal to the 
real dimension of V at x. 

Note that in the above definition, we consider a polynomial with coefficients in Q(C)- We state 
below that for a generic choice of an n x n invertible matrix A with entries in Q, / satisfies N. 

Indeed we need such a statement because our algorithm performs symbolic manipulations on 
the input by introducing an infinitesimal Q (to reduce the study to bounded semi- algebraic sets) 
and next chooses randomly A to ensure that after applying the change of variables x — t- A~^x 
some polynomial satisfies N. 

Thus, the rest of this Section is devoted to prove that up to a generic choice of A G GL„(Q), 
f^ satisfies N. 

Proposition 2. There exists a non-empty Zariski open set G C GL„(C) such that for A G (^ n 
GL„(Q), /^ satisfies N. 

Proof. We will prove that there exists a non-empty Zariski open set ff' C GL„(C) (resp. ff" C 
GL„(C)) such that for A G ^' n GL„(Q), /^ satisfies Ni (resp. N2). Taking ff = ff' n G" is 
sufficient to conclude. 

We start with Ni. [3, Proposition 3] (see also [5, Theorem 6] for a more general statement) 
states that when / has coefficients in Q and defines a smooth algebraic set ^ C C", there exists 
a non-empty Zariski open set il' C GLn(C) such that for all A G GLn(Q) n 0' and 1 < i < ra, 
crit(7rj, V^^ is either empty or is equidimensional of Krull dimension i — \. The proof of this result 
is based on the use of the Weak Transversality Theorem of Thom-Sard (see e.g. [27]). It allows 
to characterize the set of "bad" matrices, i.e. the complement of ^ in GL„(C) as the smallest 
algebraic set containing the critical values of a polynomial mapping with coefficients in the same 
base field as the one containing the coefficients of /. 

By [42, Lemma 3.5], 14 is smooth. Thus, one can apply mutatis mutandis the proof of [3, 
Proposition 3] to f — £ with C(C)(e) as a base field. We obtain the existence of a non-empty Zariski 
open set O' C GL„(C(C)(e)) such that for A G fi' and 1 < i < n, CTit{7ri,Vj^) is either empty 
or is equidimensional of Krull dimension i — 1. Recall that GL„(C(C)(e)) — ^' is Zariski closed 
and is characterized as the set of critical values of a polynomial mapping in Q{£,Q since f — £ 
has coefficients in Q(e, (^). If we multiply this polynomial by the least common multiple of the 
denominators of its coefficients, we obtain a polynomial P with coefficients in Q[e, Q. Without loss 
of generality, we can assume that the coefficients of P have no non-trivial gcd. Hence, P can be 
written as 

P = Po + £''^Qe + C^Q( + eCQ , 

where v^ and V(^ are positive integers and 



• Pq has coefficients in Q (it is obtained by substituting e and C by in P); 

• Qe (resp. Q(^) has coefficients in Q[e] (resp. Q[C]) and is not identically when e (resp. C) is 
substituted to 0; 

• Q has coefficients in Q[e,C]. 

Note that since the coefficients of P have no non-trivial gcd, at least one of the three polynomials 
Pq, Qe and (5( are not identically 0. If Pq 7^ (resp. Q^ 7^ or Q^ 7^ 0), we define 0^' C GL.„(C) as 
the non-empty Zariski defined by Pq 7^ (resp. Q^ 7^ or Q^ 7^ 0). 

Now remark that since £ and C, are infinitesimals, for all 0'' n GL„(Q) C il'; in other words for 
all A G ^' n GL„(Q), /^ satisfies Ni as requested. 

Now we deal with N2. Below, by abuse of notations the extensions of cartesian products ]0, 1[* 
in R(C) are denoted by ]0, 1[; also by convention ]0, 1[*^= {0}. The set V {^ R(C)" is semi- algebraic 
and so we can partition it into smooth semi-algebraic sets Si,..., Si [13, Chapter 5, Section 5], 
and homeomorphic to ]0, 1[^, . . . ,]0, Ip, respectively, where < A; < n. By [18, pp. 47], for 
X G ypi R(C)"5 the local real dimension of F f] R(C)" at x is given by max^^^dj, where Si denotes 
the Euclidean closure of Si . 

For 1 < i < £, we denote by Vi the Zariski closure oi Si. By [19, Proposition 2.8.2] note that the 
KruU dimension of Vi \s di, for 1 < i < i. By Noether normalization [2], there exists a non-empty 
Zariski open set Of C GL„(C(C)(e)) such that for all A e Vi'l and x G C(C)'^S 7rj/(x) n Vf' is 

finite. As a consequence, for all x e R(C)*, ^d/(x) n Si is finite. We let 9." = flLi^i'- The 
complement of Vt" in GLn(C((^)(e)) can be characterized using algebraic elimination algorithms 
(such as Grobner bases) run with parameters as for entries of a generic matrix; we refer to [37] for 
such algebraic algorithms for detecting defects of Noether normalization. As above, since the input 
has coefficients in Q(C)) 0" can be defined by a polynomial inequality P 7^ with coefficients in 
Q[e, C]: Following mutatis mutandis the approach for proving Ni, one can deduce from the inequality 
P 7^ another inequality with coefficients in Q defining a non-empty Zariski open set ff" such that 
0" n GL„(Q) C 9.". This finishes the proof. D 

4 Geometric Statement 

In this Section we let / and 51, . . . , ^s be polynomials in Z[Ai, . . . , A„], 5 C R" be the semi- algebraic 
set defined by 

/ = 0, 51 >0,...,5s >o 

and [/ C R" be the open semi-algebraic set defined by 51 > 0, ... ,5^ > 0. We will also consider 
the algebraic set F C C" defined by / = 0. Given an infinitesimal e, we denote by Ve C C(e)" the 
algebraic set defined by / = e. 

The rest of this Section is devoted to prove this result below. 

Proposition 3. Assum,e that f satisfies N, is non-negative over R", and that V^ fl R" is bounded. 
Then for < i < n, (lim£_j.o crit(7rj, Vs)) f]U = S if and only if the real dimension of S is < i — 1. 

Remark 4. [6, Proposition 3.4] provides a similar statement but with different assumptions on f 
that are not suitable for our setting. 

Before proving the above result, we start with a few lemmata. 



4.1 Auxiliary results 

The following proposition is a variant of statements that are commonly used in computational real 
algebraic geometry (see e.g. [15, 42] or [13, Proposition 12.38]). 
We let C be an infinitesimal. 

Proposition 5. Assume that f is non-negative over R" and that y n R" is non-empty. Let C be 
a semi- algebraically connected component of V CiR^. Then there exist semi- algebraically connected 
components C^^i, . . • , Cs/ of Vsf] R(e)"', such that C = IJj=i li™e-5>o Cs,i- 

Moreover if there exists a ball i? C R" such that C <Z B and C does not intersect the boundary of 
B, then C^^i C ext(-B, R(e)) and C^^i does not intersect the boundary of ext{B, R(e)), for 1 < i < i. 

Proof. Consider x G C. By assumption, / is non- negative over R". Since /(x) = there exists a 
semi- algebraically connected component 5 C R" of the semi-algebraic set defined by / > such 
that X is in the closure of S (for the Euclidean topology). Then, for all r > the ball B{x,r) 
centered at x of radius r contains a point of S (at which / is positive). By the curve selection 
Lemma [13, Theorem 3.19] there exists a continuous semi-algebraic function 7 : [0, 1] — t- 5 with 
7(0) = X and /(7(i)) > for t / 0. 

Consider the extensions 7 = ext(7, R(e)) and / = ext(/, R(e)). By the semi-algebraic interme- 
diate value Theorem [13, Theorem 2.11], there exists tg S ext([0, 1], R(e)) such that /(7(te)) = e. 
We denote ^{te) by x^ and we have lim£_^oXe = x because x^ G ext(i?(x, r), R(e)) for all r > 0. 
Also, let Cxj be the semi- algebraically connected component of 14 n R(e)'^ that contains x^; we 
associate to x this semi-algebraically connected component Cx^- of 14 n R(e)". 

Since there are finitely many semi-algebraically connected components of 14 n R(e)", there are 
finitely many semi-algebraically connected components of Cs,i, . . . , Cs/ which are associated to a 
point X in C. This proves that C C Uj=i li™£-s>o Ce,i. 

We prove now that IJi=i l™£-s>0 C'e.i C C. Let x G lim^-^o Cs,i for some 1 < i < £; we need 
to prove that x G C. Then, there exists Xg G Cs,i bounded over R such that lim£_!.o x^ = x. 
Now, recall that Cg^j is associated to some point x' G C; this means that there exists x^ G Cs^i 
bounded over R such that lime_j.ox^ = x'. Since Ci^s is semi-algebraically connected, there exists 
a continuous semi-algebraic function 7 : [0, 1] — )• Cj^^ such that 7(0) = x^ and 7(1) = x^; we have 
that r = 7([0, 1]) is a semi-algebraically connected semi-algebraic set. By [13, Theorem 3.20], 
r = 7([0, 1]) is closed and bounded. By [13, Proposition 12.36], we conclude that lim£_j.or is semi- 
algebraically connected. Since lim£_^o is a ring homomorphism, lim£_^or is contained in 1/ n R"". 
Now notice that lim£_!.o(7(0)) = x and that lim£^o(7(l)) = x'. Since we have proved that lim£_!.o T 
is semi-algebraically connected we deduce that x G C as requested. 

Now, we assume that C is bounded and let i? C R" be a ball such that C d B and C does 
not intersect the boundary of B. To conclude the proof it remains to prove that for 1 < i < ^, 
Cs,i C ext(i?, R(e)) and Ce^i does not intersect the boundary of ext(i?, R(e)). 

Consider x^ G Ce^i bounded over R such that lim£_>.o x^ G C. Such a point exists, as we argued 
in the first part of the proof that lime_^o C'e.i C C. 

We argue by contradiction. Assume that there exists x^ in the Euclidean closure of C^^i — 
ext(i?, R(e)). Since Cs^i is semi-algebraically connected, there exists a continuous semi- algebraic 
function 7 : [0, 1] — )■ Ce,i such that 7(0) = x^ and 7(1) = x'^. Note that by the intermediate value 
theorem [13, Theorem 2.11] (applied to the polynomial defining the boundary of B) there exists 
"d G [0, 1] such that 7(1?) lies in the boundary of ext(S, R(e)). We deduce that lim^-^o 7(''?) belongs 
to the boundary of B. 
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By [13, Theorem 3.20], 7([0, 1]) is closed and bounded. We also notice that 7([0,1]) is semi- 
algebraically connected. By [13, Proposition 12.36], we conclude that lime_!.o7([0, 1]) is semi- 
algebraically connected. We deduce that lim£_i.o(7([0, 1])) C C. We deduce that lim£_j.o7(i?) 
belongs to C. 

Consequently, we have Hui^^q^^Il)) belongs to C and to the boundary of B. This contradicts 
the fact that, by assumption, C does not intersect the boundary of B. Thus, we conclude that 
Ce^i C ext{B, R(e)) and C^^i does not intersect the boundary of ext(i?, R(e)). □ 

The following lemma relates the real dimension of a semi-algebraic set in R(e)" with the real 
dimension of its image by lim£_j>o- Its proof uses quite standard tools from real algebraic geometry. 

Lemma 6. Let S^ C R(e)"' (resp. Z^ C C(e)"'j he a semi- algebraic (resp. constructible) set and 
d he its real (resp. Krull) dimension. Then lime_>o Se (resp. Ivch^^q Z^) has real (resp. Krull) 
dimension less or equal to d. 

Proof. The proof of this statement relies on [13, Proposition 5.29] which states that if A is a semi- 
algebraic subset of R™ and /i : ^ — )• R" is a semi-algebraic mapping, then the real dimension of 
h{A) is less than or equal to the real dimension of A. Thus, it suffices to prove that the ring 
homomorphism lime_s>o is a semi- algebraic function, i.e. a function whose graph is a semi- algebraic 
function. This is a quite routine statement in real algebraic geometry; we give the proof since 
we could not find a reference stating it explicitly. To do that, we reuse some ingredients of [13, 
Proposition 12.36]. 

Recall that by assumption Se is a semi- algebraic set of R(e)'^ and that, by definition R(e) is the 
real closure of R(e). Then, by [13, Proposition 2.82], there exists a quantifier- free Boolean formula 
of conjunctions and disjunctions of polynomials in R[e][Xi, . . . ,X„] which define S^- Below, we 
denote by ^(X, e) such a formula. 

Consider the following set 

T = {(x, Xn+i) G R"+^ I l-Cx, Xn+l) A 2;„+i > 0} , 

and let T its closure and H be the hyperplane defined by Xn+i = in R"'^"^. As in the proof of 
[13, Proposition 12.36], the following equalities hold: 

S = \\m.Se = Tr\H. 

We can express the limit using a formula in the first order theory of reals. For x G iSg and y G 5 
we have 

$1 := [(Vr > 0) (3eo) : (Ve) (0 < e < eo) ^ ||x - y f < r^] 4^ lim x = y . 

By quantifier elimination over the reals (see [24, 49] or [13, Theorem 2.77]), there exists a quantifier 
free formula ^i which is equivalent to 'I'l. The set 

{(x,y)GR(e)"xR"|vI/A(yGcS)AvI/i} , 

is the graph of lim£_>o ^md is semi- algebraic. This concludes the proof. 

The proof in the complex case uses exactly the same techniques as above transposed to alge- 
braically closed fields, i.e. if A and B are constructible sets and h : A ^ B \s a regular map (see 
[46, Chapter 1]), then the Krull dimension of h{A) is less than or equal to h[B) and quantifier 
elimination over algebraically closed fields (see e.g. [13, Chapter 1]). D 
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Lemma 7. Consider x G Fn R" and assume that f is non-negative over R" and that there exists a 
neighbourhood ^ C R" o/x such that y P| =^ P| R" zs a finite set. Then, there exists x^ G crit(7ri, V^) 
such that lim£_j>o x^ = x. 

Proof. By assumption yp|=^P|R" is a finite set. Hence, there exists a sufficiently small r E R, 
r > 0, such that / is positive at all points in a ball, B, with center at x and radius r, except x. 
That is X is the only point of V{^R^ in B] hence {x} is a bounded semi-algebraically connected 
component of F P| R". 

By Proposition 5, there exist semi-algebraically connected components Ce,i, . . . , Ce/ of 14, R(e)" 
such that {x} = Ui=i lhxL£_^o C'e,*, Ce,i C ext(S, R(e)), and Ce,i does not intersect the boundary of 
ext(5, R(e)), 1 <i <^. 

We deduce that for 1 < i < ^, Ce^i is closed and bounded. However, Cg/i C ext(i?, R(e)) 
and it is closed and bounded; as a consequence, it has a non-empty intersection with crit(7ri,^). 
Since we already observed that lime_s.o Ce^j C {x}, we deduce that for all x^ G ci\i{'Ki.,Vir) n Ce^j, 
limg^o Xg = X. Our conclusion follows. D 

Lemma 8. Assume that f is non-negative over R" and that y n R" is bounded and non-empty. 
Then, for 1 < i < n, crit(7rj, Vs) is not empty and intersects all bounded semi-algebraically connected 
components ofVeO R(e)". 

Proof. Since V OR^ is bounded and non-empty there exists a semi-algebraically connected com- 
ponent C and a ball S C R" such that C C B and C Ci B = 0. By Proposition 5, we deduce 
that there exist semi-algebraically connected components C^^i, . . . , Ce^fc of 14 n R(e)" such that 
C = U£=i l™e->-o Ce/. Moreover since we assume that ^ n R" is bounded, there exists a ball B C R^ 
such that C C B and C Ci B = ^. Using again Proposition 5, we conclude that C^^^ C ext(5, R{e)) 
for 1 < £ < k; hence Cs/ is closed and bounded. 

Below, we prove that any bounded semi-algebraically connected component C^ of V^nR(e)" has 
a non-empty intersection with crit(7ri, 14). Then 7ri(Ce) is closed and bounded (see [13, Theorem 
3.20]). The extreme values of 7ri(C£) are attained at critical points of the restriction of tti to 
Vs n R(e)". In other words, C^ has a non-empty intersection with crit(7ri, 14). Now, note that 
crit(7ri, 14) C crit(7rj, I4), by definition. D 

Lemma 9. Let f G C[Xi, . . . , X„] and V C C^ be defined by f = 0. Take a = (ai, . . . , aj„i) € C* 
and let Vi^a the algebraic set V f~\7r^_i{a) and cpi be the canonical projection (xi,...,x„) 1— t- Xj. 
Then, the following holds: 

cnt{ipi,Vi^a) C crit(7rj,y). 

Proof. By definition crit(7rj,l/) is defined by the vanishing of / and of the partial derivatives 
cDc~' • ■ ■ ' dx- Besides, the set Vi^a is defined by the system 

f = 0,Xi-ai = ■■■ = Xi^i - aj_i = 0. 

By definition, crit(c/3i, 1^^^) is defined by the above equations and the vanishing of the maximal 
minors of the Jacobian matrix associated to f,Xi — ai, . . . , Xj_i — aj_i,Xj. The triangular shape of 
this Jacobian matrix implies that cnt{ipi, l^,a) is defined by / = 0, Xi —ai = --- = Xj_i — ai-i = 
and the vanishing of the partial derivatives q-^^ , • • • , qx ■ 

The lemma follows from the observation that the system that defines crit(7rj, Vi) is contained in 
the system that defines ciit{ipi,Vi^a)- D 
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4.2 Proof of Proposition 3 and Consequences 

Proof of Proposition 3. First we prove the necessary condition. Assume that for a given i E 
{0, . . . , n} it holds 

(lim(crit(7ri,y,)))f|C/ = 5 . 

We prove that this imphes that the real dimension of S" is < i — 1. When S is empty, this is 
immediate; thus we can assume that S is non-empty. 

By definition, V f]U = S, therefore 1^ f] R*^ is non-empty; note also that S is bounded since we 
assume that F Q R*^ is bounded. 

Combining Lemma 8 and Ni, we conclude that crit(7rj, Vs) has Krull dimension i — 1. Lemma 6 
implies that lim£^o(ci'it('/rj, Ve)) has Krull dimension < i—1. We deduce that (lim£^o(crit(7ri, 14))) P| R" 
has real dimension < i — 1. Consequently, (lime_i.o(crit(7rj, V^))) P| L'^ has real dimension < i — 1. 
By assumption, (lime_s.o(crit(7rj, Vs))) f]U = S; we conclude that S has real dimension < i — 1. 

Next, we assume that the real dimension of 5" is < i — 1; we prove below that this implies that 

(lim,^o(crit(7ri,K)))nf^ = '5- 

The inclusion (lime_s.o(crit(7rj, T^))) P| [/ C 5* follows from the fact that crit(7rj,F£) C Vg, 

liuie^o VeCV,andVf]U C S. 

It remains to prove the inverse inclusion, that is 5 C (lims^Q {cr it {it i, V^))) f] U. 

If S is empty, then it is immediate that S C (lim£_^o(crit(7rj, V^))) f] U. In the sequel, we assume 
that S is not empty. 

Let X = (qi, . . . ,a„) G S; since S = U f]V hy definition, we have x G [/ and x £ V f] R". 
Thus, we need to prove that there exists a point x^ G crit(7ri, Vs), such that liuis^o x^ = x. This is 
what we do below. 

Recall that we have assumed that the real dimension of S" is < i — 1. Since V f]U = S, we 
deduce that the real dimension of y p| [/ is < i — 1 . It follows that the local real dimension of 
y Pi [/ at X is at most i — 1. We claim that the real dimension of y f] R" at x is < i — 1. 

Indeed, let Si,. . . ,Sk be smooth semi- algebraic sets of respective dimension dj which form a 
finite partition of V. The local real dimension of y P| R" at x is the maximum of the dj's for j 
such that X G Sj. Since U is open, Sj f]U is either empty or smooth of dimension dj. Note also 
that the semi-algebraic sets Sif]U,...,SkC\U form a partition ofVf]U. Therefore, the local real 
dimension of ^ f] R" at x is the same as the local real dimension ol V f]U at x which is < i — 1 as 
requested. 

Denoting 7rj_i(x) by Xj_i, we deduce by N2 that -K~_-^^{-Xi-i)f]V (^R"" is finite. Now, let / be 
the polynomial obtained after instantiating the first i — 1 variables with the first i — 1 coordinates 
of X, that is 

/ = f{Xi,. .. ,Xn) = f{ai,... ,ai-i,Xi,. . . ,Xn) ■ 

We let X = {oi, . . . , On) and V (resp. Vs C C(e)"~*+-'^) be the algebraic set defined by / = 
(resp. f = e). By assumption, / is non- negative over R". Therefore, / is non- negative over R""***"^. 
We also consider the canonical projections 

(^j : (xi,...,x„) -^ Xj and y?i : (xi, . . . ,x„) -^ Xj. 

By applying Lemma 7 to F and /, there exists Xg G crit((^j, Vs), such that lim£_s.o Xg = (a,, . . . , an)- 
Now define x^ = (ai, . . . , ai_i, x^) and V^ = T^ f| 7rr_\(ai, . . . , ai_i). 
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Since x^ G crit((^j, 14), it is immediate that x^ G cnt{ipi,V^), and using Lemma 9 x^ S 
cvit{7ri,Vs). To summarize, we have estabhshed x^ G crit(7ri,V^) and limg-j^oXj = x as requested. 
This finishes the proof. D 

The following corollary is a direct consequence of Proposition 3 and it is crucial for the proof 
of correctness of our algorithm. Recall that, by convention, crit(7ro, V^) is the empty set. 

Corollary 10. Assume that f satisfies property N and that V is bounded. Let d he the real 
dimension of S. Then, d > if and only if there exist an integer i in {l,...,n} such that 
(limj_!.o crit(7rj, Vs)) P\U ^ (lim^-^o crit(7ri_i, V^)) f] U and d + 1 is the largest of these integers. 

Proof. Note that the assumptions of Proposition 3 are satisfied. This implies that, for < z < n, 
(lim^^o crit(7ri, Ve))f]U = Sif and only if d < i — 1. As a consequence, we deduce that 

S = (lim crit(^d+i, F,)) Pi [/=... = (lim crit(^„, F,)) f| U 
and S / (lime_>o crit(7rj, 14)) P| [/ for 1 < i < d. This concludes the proof. D 

5 Algorithm and proof of Theorem 1 

5.1 Main algorithm 

The input of our algorithm is two sets of polynomials F = (/i, . . . , fp) and G = (51, . . . ,5s) in 
Z[Xi, . . . , Xn]. The output is the real dimension of the semi- algebraic set ^ C R" defined by 

/i = • • • = /p = 0, gi>0,...,gs>0 . 

The main idea of the algorithm is to perform symbolic manipulations to define a polynomial / 
with coefficients in R[C] such that 

• / and the semi- algebraic set S C R(C)"^^ defined by / = and gi > 0, . . . ,gs > satisfy the 
assumptions of Corollary 10, and 

• the real dimension of S is the same as the real dimension of J/'. 

Then, we apply the results of the previous section (Section 4) using as ground field (that is the 
field where the coefficients of the input polynomials belong to) R(C)- 

Let us recall the notations that we use: S C R(C)"^^ is the semi-algebraic set defined by 
f = 0, gi > 0, . . . , gs > 0, U C R(C)"^^ is the open semi- algebraic set defined by gi > 0, . . . , ^f^ > 0, 
and V C C{C)"'~^^, resp. Ve C C{C) {e)'^~^^ , is the algebraic set defined by / = 0, resp. f = e. 

Prom Corollary 10, d > if and only if there exists an integer i in {1, . . . , n} such that 



(lim crit(^i, Ve)) fl [/ / (lim crit(^i_i, Ve)) f| U 



and d+1 is the largest of these integers (notice that by convention crit(7ro, K) = 0; see Section 2). 
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Subroutines. We need three subroutines. Tlie first one is called Random; it takes as input an 
integer n and returns a randomly chosen n x n matrix in GL„(Z). 

The second routine IsEmpty takes as input a polynomial H and a set of polynomials Q = 
(Gi, . . . ,Gs) in Z[(,r],e][Xi, . . . ,Xn], where C, V^ ^i^d ^ ^^^ infinitesimals. Let V C C(C)(??)(e)" be 
the algebraic set defined hy H = and U G R{C) [r]) (e)^ the open semi-algebraic set defined by 
Gi > 0, . . . , Gs > 0. The subroutine decides if the semi- algebraic set V f]U is empty or not. 

Let 5 be the maximum of the degrees of the monomials in ^,rj,e,Xi, . . . ,X„ appearing in H 
and Q. If Z = Z, then we denote by r the maximum bit size of the integers appearing in H and 
Q. The algorithm returns True if the semi-algebraic set V f^U is empty, and False otherwise. It is 
based on [13, Algorithm 13.1]. 

Lemma 11. Using the above notations, algorithm, \sErr\pty(H,Q) decides if the semi-algebraic set 
defined by H = and Gi > 0, . . . ,Gs > is empty or not within {s 5)^^' arithmetic operations in 
Z. When Z = Z, then the Boolean complexity o/ IsEmpty is r(s(5) '"^. 

When IsEmpty is called by the main algorithm, the input polynomials have coefficients in Z[Q. 

The third subroutine, DisjointPolar, takes as input a polynomial / G Z[C][Xi, . . . ,Xn+i], G = 
(gi, ■■■,9s) C Z[Xi, . . .,Xn], Q G Z[C][Xi, . . . ,X„+i], and an integer i G {1, . . . ,n} (where C is an 
infinitesimal that we manipulate as a variable). The polynomial Q defines a sphere in R((^)"'"'"^ that 
strictly contains the semi- algebraic set S. For example Q = {C{Xf -|- • • • -|- X^ + ^n+i) ~ 2)- 

The routine DisjointPolar(/, G, Q,i) returns True if 



flim crit(7ri, Ve) n C/ W f lim crit(7ri_i, Ve)) n U, 



and False otherwise. Below, we denote by 6 the maximal degree of the monomials in (, Xi, . . . , Xn 
appearing in /, G and Q. 

The routine DisjointPolar is described Section 5.2, where we also prove the following Lemma. 

Lemma 12. Let f e Z[C][Xi, . . . ,X„+i], G = (51, ■.■,9s) C Z[Xi, . . . ,X„], B G Z[C][Xi, . . . ,X„+i], 
and i > 0. Let S be the semi- algebraic set defined by {f = A G > 0}, and let Q < define a 
ball that strictly contains S. Assume that f satisfies N and is non-negative over R(C)" O'lT-d that 
(lim£_s.ocrit(7rj, 14)) Pi [/ is not empty. There exists an algorithm DisjointPolar which decides if 
(limj_!.o crit(7ri, V^)) nU j^ (lime_>o crit(7rj_i, V^)) D U within {s 5)*^(") operations in Z. 

// Z = Z and the maximum bit size of the integers in the coefficients of the input polynomials 
is T, then the Boolean complexity of the algorithm is r {s5)^^'^' . 
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We can now describe our main algorithm ComputeRealDimension. 



Algorithm 1: ComputeRealDimension(F, G) 

Input: F = (/i, ...Jp) and G = (gi, ... ,5^) in Z[Xi, . . . ,X„] 

Output: The real dimension of the semi-algebraic set defined by /i = • • • = /p = 0, gi > 0, . . . , g^ > 

2 if lsEmpty(/o, G) then return —1 ; 

3 / ^ /o + iCiXf + • • ■ + X2 + X2^i) - 1)2 ; 

4 Q^(CiX! + --- + Xl + X-^+,)-2); 

5 A <— Random(n) ; 

6 /^(X) ^ /(AX) and G^(X) 4- G(AX) and Q^(X) ^ Q(AX); 

7 for 71 > j > f do 

8 1^ if DisjointPolar(/^, G^, Q^, i) then return i — I ; 

The proof of Theorem 1 consists of proving the correctness and complexity estimate of Com- 
puteRealDimension. 

Correctness. Let d be the real dimension of the semi-algebraic set =5^ C R" defined by /i = 
■■■ = fp = 0,gi>0,...,gs>0. 

At Step 2 we test whether the semi- algebraic set is empty or not. In the sequel we assume that 
the semi- algebraic set is not empty, and so d > 0. 

We also denote by ^ C C(C)""^^ the algebraic set defined by / = where / is defined at 
Step 3. Assume for the moment that the semi-algebraic set S has real dimension d and that 
V n R(C)" is bounded. At Step 6, the real dimension of the semi-algebraic set S-^, defined by 
/ = and §1 > 0, . . . , g^ > 0, is also d. Since A is randomly chosen we can assume that it lies 
in the non-empty Zariski open set iff defined in Proposition 2. Therefore, / satisfies property N 
(see Definition 2) and V^ n R(C)"+^ (and consequently S-^) is bounded. In other words, all the 
assumptions of Corollary 10 are satisfied. 

In the for-loop, starting with i = n, the algorithm looks for the largest integer i such that 

f lim crit(7ri, Vj^)) nU^ ^ f lim crit(7ri_i, Vj^)) H U^ . 
\e—>-0 ) \e— >0 / 

Each time we enter in the loop (lime_!.o crit(7rj, l/^^^)) n \J^ is not empty (and actually equals 
S'^). The other assumptions of Lemma 12 are obviously satisfied. Also, by Corollary 10 and 
Lemma 12, we deduce that when d > 0, at Step 8 the algorithm will return d. 

To finish the proof it remains to establish that S is bounded and its real dimension is the same 
as the real dimension of 5^ . 

We start with the boundedness statement. The polynomial / is the sum of the squares of 
fi, . . . , fp and the square of the polynomial C(-'^i + • • • + X^ + ^n+i) ~ ^- The set of roots in 
R(C)"'"'"^ of this latter polynomial is the n-dimensional sphere, Sn+i) with center at the origin and 
radius l/i/C- It is straightforward that V n R(C)"+^ is bounded and since S = V CiU, we deduce 
that S is bounded. 

Now, we prove that S has the same real dimension, d, as the semi- algebraic set J^. The proof 
consists of two steps. We prove consecutively that 
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(i) ext(^, R(C)) C R(C)" has real dimension d, and 

(ii) S has the same real dimension as ext(=5^, R(C)) C R(C)"'- 

Statement (i) is a straightforward consequence of the Lemma below. 

Lemma 13. Let S be a semi-algebraic subset of R" and consider its extension ext(5, R(C))- Then 
S and ext(5, R(C)") have the same real dimension. 

Proof. Let d be the real dimension of S and d' be the real dimension of ext(5, R(C))- Then there 
is an injective map (/J : 5 — )• [0, 1]"^. We consider the extension of ip, that is ext((^, R(C))- It is also 
injective [13, Exercise 2.17] and by the definition of the real dimension (see Definition 1) we deduce 
that d' < d. Now remark that lim^_!.oext(5, R(C)) = S. By Lemma 6 we deduce that d < d' which 
finishes the proof. D 

To prove (ii), we consider the restriction of 7r„ to S. It is a semi-algebraic and injective function, 
since for all (x, x^+i) G R(C)" ^ R(C) with x G exi{.y,R{()) and (x, x„_t-i) G Sn+i> Xn+i = 
Jj-||x||2. Moreover, 7r„(5) = ext(^, R(C)). 

By [13, Proposition 5.29] we deduce that the real dimension of S is equal to the real dimension 
of ext(=y, R(C)); which is d by (i). D 

Complexity Analysis. The complexity estimate is straightforward by Lemmata 11 and 12. The 
first call to IsEmpty costs at most [s D)^^^' arithmetic operations. Next, the algorithm calls n times 
the sub-routine DisjointPolar with input polynomials of total degree 0{D). Each call costs {s D)^"^' 
operations in Z. Since all functions in the complexity class n{sD)'^^^> lie in the complexity class 
{s D)^""' , we deduce that ComputeRealDimension runs within {s D)^^' arithmetic operations in Z. 
This is a probabilistic bound because of the random change of coordinates that we apply to / 
and G at Step 6 to ensure property N (Definition 2). 

5.2 The subroutine DisjointPolar and Proof of Lemma 12 

In order to describe DisjointPolar, we need to recall some fundamental algorithmic specifications 
and complexity results in computational real algebraic geometry. Sometimes we need to prove some 
statements which are quite folklore in the area; proofs of these facts are given in the Appendix. 

Quantifier Elimination over the reals. Let X = (Ai, . . . , A„) and Y = (Yi, . . . , Vj) be finite 
sequences of variables, Jl be the existential quantifier 3 or the universal quantifier V and 'P(X, Y) 
be a Boolean function of s atomic predicates ifj(X,Y) Oj 0, where Oj G {>,<,= } and 
Hi G Z[X,Y], for 1 < i < s. We let 6 = max(deg(-fri), 1 < i < s) and when Z = Z, r is the 
maximum of the bit size of the coefficients in the Hi 's 

One-block quantifier elimination over the reals consists in computing a quantifier-free formula 
which is equivalent to the first order quantified formula ^ : (OX G R")P(X, Y). 

The following theorem is a simplification of the general purpose quantifier elimination algorithm 
in [13, Algorithm 14.5 and Theorem 14.16]. 

Theorem 14 (One Block Quantifier Elimination over the reals). There exists an algorithm 
OneBlockQuantifierElimination which takes as input ri,P,X, Y and returns a quantifier free formula 
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^ of the form Vi=i Aj=i hij >ij 0, where hij G Z[Y'i, . . . ,Yt] and \>ij G {>, =}, that is equivalent 
to <I>, such that 

/f Z = Z, i/ien t/ie maximum bit size of the coefficients of hij is bounded by t6^^^^'^"'' . The above 
transformation requires s(*+i)("+i)(^(*+i)0(n) arithmetic operations in Z. The Boolean complexity of 
the algorithm is rs(*+i)("+i)5(*+i)^("). 

Limits of semi-algebraic sets in R{0{r]){e) when e — )• 0. The routine UnivariateLimit takes 
as input a quantifier-free formula ^, which is a disjunction of conjunctions of polynomial equa- 
tions/inequalities in Z[C,r/,e][Z], and the infinitesimal e. This formula defines a semi- algebraic set 
5^ in R{(){r]){£). It outputs a semi-algebraic description of lime_5.o5e. This is based on quantifier 
elimination over the reals [13, Theorem 14.16]. 

Lemma 15. Let 5 be the maximum of the degrees of the momonials in CiV^^ of the polynomials in 
^ and i be the number of polynomials in ^. One can compute UnivariateLimit(^,e) within [£S)^^' 
arithmetic operations in Z. 

When Z = Z and t is a bound on the bit size of the integers of the coefficients in ^, then 
UnivariateLimit(^,e) runs within T(id)^^^' bit operations. 

Computing limits of critical points. The routine LimitsOfCriticalPoints takes as input a 
polynomial H and a set of polynomials G = {Gi, . . . ,Gs) in Z[C][Xi, . . . ,X„]. We denote by 
Vs C C(C)(e)" the algebraic set defined by i:f — e = 0, by C/ C R(C)" the open semi- algebraic 
set defined by Gi > 0, . . . ,Gs > and by 6 the maximum of the degrees of the monomials in 
C, Xi, . . . , Xn appearing in H and Q. If Z = Z, then we denote by r the maximum bit size of the 
integers appearing in H and Q. It is based on [13, Algorithm 13.1] and [13, Algorithm 11.20]. 

Lemma 16. Assume that H is non-negative over R(C)" o^nd that crit(7ri, 14) is finite. There exists 
an algorithm LimitsOfOfCriticalPoints which takes as input H and Q as above, and returns True if 
(lim£_j.o crit(7ri, 14)) Pi [/ is non-empty else it returns False, within {sS)^^"^' arithmetic operations 
in Z. 

When Z = Z and r is the maximum bit size of the coefficients, then the Boolean complexity of 
LimitsOfCriticalPoints is t{s5)^'^"-\ 

The routine IsRealizable. A typical call of this routine is lsRealizable($), where $ is a union of 
disjunctions of univariate polynomials in Z with coefficients in Z[(^, 7;]. 

The subroutine calls Algorithm 10.13 from [13] to compute all realizable sign condition of the 
polynomials and checks if there is at least one sign condition that is compatible with <I>. In this 
case it returns True, otherwise it returns False. 

Lemma 17. Let ^ be a union of disjunctions of univariate polynomials in Z with coefficients in 
Z[C, r/] such that all monomials in Z,C,.,rj appearing in $ have degree at most 5. 

The complexity 0/ lsRealizable($) is 16^^^' operations in Z. If 7. = Z and the maximum bit size 
of the coefficients of the polynomials of the input is t, then the Boolean complexity is t 1 5^^^' . 
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The routine DisjointPolar. Let / G Z[Q[Xi, . . . , Xn], a set of polynomials G = {gi,...,gs) C 
Z[Xi, . . . , Xn], a polynomial Q G Z[(] [Xi, . . . , Xn] and an integer i, such that 1 < i < n. 

The open semi- algebraic set of R(C)" defined by (71 > 0, . . . ,5^ > is denoted by U. We 
introduce another infinitesimal r], with (^ > 77 > 0, and we consider the semi- algebraic set Urj C 
RiOif])"" defined hy gi > r], ... ,gs > r]. We also denote the latter inequalities hy G > rj. Note that 
[/^ is closed. 

In the sequel, we denote R(C)(ry) by R and C((")(^) by C. 

Finally, we introduce one more infinitesimal e with Q > -q > e. Let V^ C C(e)" be the algebraic 
variety defined by / — e = 0. To make the notation simpler, ext(C/^,R(e)) will be denoted by Un,e 
and ext(i?,R(e)) will be denoted by S^. 

We assume that 

• / satisfies property N (Definition 2) and is non-negative over R(C)"'; 

• that Q < defines a ball B that strictly contains the semi- algebraic set defined by / = and 
51 > 0, . . . , 5-5 > 0. 

• and the solution set of gi > 0, . . . , (7^ > has a non-empty intersection with crit(7rj, V^); 

Then, DisjointPolar(/, G, Q, i) returns True if (limg^o crit(7rj, Ve)) fl ^ / (lim£^o crit(7rj_i, V^)) H C/, 
and False otherwise. Before giving a detailed description, we briefiy give a geometric view of the 
operations performed by DisjointPolar. 

The algorithm works as follows. We consider separately the case i = 1 (Step I). By convention 
crit(7ro, Ve) is the empty set. Therefore, (lim£_^o crit(7ro, Ve)) {^U = % and if i = 1, then it suffices 
to check whether (lime_j.o crit(7ri, Ve)) {^U is empty or not. 

When i > 2, the algorithm starts by testing if crit(7rj_i, 14) P| C/^^e is empty. We will see that 
when this is the case, it implies that (linie^o crit(7rj_i, V^)) P| [/ is empty. Since by assumption 
(lime_^o crit(7rj_i, Ve)){^U is not empty, the routine returns True. 

When crit(7rj_i, Ve) H C/r?,e is empty, the routine constructs a formula that defines the following 
semi- algebraic set 

yl, = {(x,y,z)ER(e)2"+i I xG(crit(7r„I4)-crit(7ri_i,K))f|[/^,,f|i?„ 

y ^ Cmt{TT^-i,Ve)[^U^,e{~\Be, \\x - y\\ = z > ^} . 

Next, we use quantifier elimination to compute a semi- algebraic description of the semi- algebraic 
set WiHe^oT^zi-^e) C R (where ttz is the projection on the Z-coordinate) . We will prove that 

flim crit(^i, Ve)) r\U^ (lim crit(^i_i, Ve)) fl C/ ^^ lim ^z(^.) fli^ > 0} ^ . 

For the correctness proof of DisjointPolar we will need the following lemma. 

Lemma 18. Assume that f is non-negative overIC' and satisfies N. //x G (lime_i.o crit(7rj, Ve)) H R", 
then there exists x^ E crit(7rj, Ve) f] R(e)"' such that lime^o ^e = x. 

Proof. Consider the polynomial H = f + Ylk=iidx)'^ ^^'^ notice that -ff(x) = 0. Then, apply 
Proposition 5 to the semi-algebraically connected component of the real solution set oi H = 
containing x. Thus, there exists x^ S R(e)" such that H{xe) — e = and lim£_s.o x^ = x. This 
implies that x^ G crit(7rj_i, 14) f] R(e)". D 
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Algorithm 2: DisjointPolar(/, G, Q,i) 



Input: / € Z[C][Xi, ...,Xn],G= (51, ... ,5s) C Z[Xi, . . . ,X„], Q G Z[C][Xi, . . . ,X„], and i 
Output: True if (lim£_^o crit(7rj, Ve)) ClU ^ (linie^o crit(7rj„i, Vs)) H U. False otherwise. 
Assumptions: / satisfies N and is non-negative over R(C)") Q < defines a ball strictly 
containing S, (lime_>o crit(7rj, Ve)) f]U ^ 9, and 1 < i < n 

1 if i = 1 then return LimitsOfCriticalPoints(/, G) 

2P,^{{fiX)-e)'+(^y + ... + (^^Y = 0AG>vAQ<0}; 



3 P2^Subs(X = Y,{(/(X)-e)2+(^)' + ---+(^)' = 0AG>r?Ag<0}) ; 

4 if lsEmpty(P2) then return False ; 

5 $^ [A(X) A P2(Y) A ^gi /O A ||X-Y||2 > Z^ A Z > 0] ; 

6 ^ ^ OneBlockQuantifierElimination(3,<l), [X, Y], [Z, C,?7,e]) ; 

7 /* \I/ = [vl'i,...,vi/^] */ 

8 $ ^ [UnivariateLimit(^i,e),... Limit(vl'/,e)] ; 

9 for 1 < fc < / do 
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if lsRealizable(^,fc A (Z > 0)) then return True ; 



11 return False 



We can now prove Lemma 12. 

Correctness. By assumption, / is non-negative over R(C)'^ and satisfies N. Thus crit(7ri,V^) 
is finite. Moreover, all the assumptions of Lemma 16 are satisfied that implies correctness when 
i = 1. 

In what follows, we assume that i > 2. 

Note that the systems Pi and P2 at Steps 2 and 3 define respectively the sets 

Cmt{TTi,Ve)f]Ur,,ef]Be and Cmt{-Ki.i,Vs) f]Ur,,ef]Be . 

At Step 4, the call to IsEmpty returns True if and only if crit(7ri_i, Ve) f] Un,e fl Be is empty. We 
claim that if this latter set is empty then (lim£_j.o crit(7rj_i, 14)) f] C/ is empty. Since by assumption 
(lim^^o crit(7rj, Ve)) {^U \s not empty, DisjointPolar runs correctly by returning True at Step 4. Now 
we prove our claim. 

Assume by contradiction that there exists x G (linig^o crit(7rj_i, V^)) P| U . Then, by assumption 
X e i? and gi(x) > 0, . . . , fl's(x) > 0. We deduce that for 1 < A; < s, we have 9fe(x) > r/. The latter 
inequality is strict because x G R(C)"' and the coefficients of g^ lies in R{(){r]). 

Since / is non-negative over R(C)", there exists x^ E crit(7rj_i, Ve)ri'R{e)^ such that lime_j.o x^ = 
X, by Lemma 18. As a consequence, for 1 < fe < s we have gki^e) ^ V (otherwise, using lime_^0; 
this would contradict x G C/^). As a consequence, we have x^ G crit(7ri_i, Ve) f] Un,e- Finally, since 
X G (lime_5.o crit(7rj_i, Ve)) H U , it is at positive distance to the boundary of B. Since lim£_^o ^e is 
X they are infinitesimally close and x^ lies in Be at positive distance to its boundary w.r.t e. We 
conclude that there exists x^ G crit(7rj_i, Ve) H t/^,^ H Be which is a contradiction. 
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In the rest of the proof we can assume now that crit(7rj_i, V^) f] [/^ f] i?^ is not empty. 
In the discussion preceding the correctness proof we considered the semi- algebraic set 

A, = {(x,y,z)GR(e)2"+i | x G (crit(7ri, K) - crit(7r,„i, K)) f| C/^,,f| S„ 

\\x — y\\ = z > 0} . 
We can also describe Af, as 

A, = {(x,y,z) G R(e)2"+Vi(x) AP2(y) A ^^^ / A ||x - y|| > z^ A z > 0} . 

We refer to Steps 2 and 3 for the definitions of Pi and P2 and the polynomials that they involve. 
If we replace the three infinitesimals e, 77, and C with three new variables e, h, and t respec- 
tively, in the formula $ defined at Step 6, then we get polynomials Pi{'X.,e,h,t), qx^ ' , and 
P2(Y, e, /i, t). In this way we define the following semi-algebraic set 

A = {(x,y,z,e,^t)GR^"+^|Pi(x,e,^t)AP2(y,e,/i,t)A ^^^^^;^'^^ / OA||x-y|| > z^Az > 0} . 
By the definition of an extension (ext), see Section 2, of a semi-algebraic set we get 

7rx,y,. (ext(^, R(e)) p|{e = e, /i = 7?, t = C}) = A, (Z R(e)" , (1) 

where 7rx,y,z : (x, y, z, e, h, t) 1— ?■ (x, y, z). We deduce that 

vr.(7rx,yz(Ext(.4, R(e)) f]{e = e,h = rj,t = (})) = ^.(A,) , 

where tTz : (x, y, z) 1— ;■ z. 

The application of OneBlockQuantifierElimination at $ (Step 6) actually computes a semi- 
algebraic formula defining TTz^(,^h,t{A). By (1) and the above discussion, one can conclude that 
substituting e, h and t by e, t] and ^ in this formula provides a semi- algebraic formula defining 

T^ziAe). 

Finally, we have to prove that indeed IsRealizable returns the correct answer. This means that 
the following formula is true: 

Bz G R, z > A z G lim TTziAs) <^ Clim crit(7ri, Vs)) f]U ^ Aim crit(7rj_i, Vs)) Q U , 

for some given value of i G {2, . . . , n}. 

The if statement. If such a z exists, then there are distinct points x^, y^ in R(e)" and z^ G R(£), 

such that (x£,y£,2;) G A^, that is 

x^ G (crit(7ri, K) - crit(7r,„i, y^)) f] U.n,e f] B^, Ye G crit(7ri_i, K) H t^r,,e fl Be , 

and lim^^o ||xe - ys\\ = z"^ > . 
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This means that x^ and y^ are not infinitesiniahy close, w.r.t. e. Notice that, by definition of A^, 
Xg and ye both lie in B^ that is bounded over R. Consequently, they are bounded over R and their 
limits as e — 7- exist. The limits as e — t- are different and thus 



lim (crit(7r„ Ve) f| U^,e f| bA / lim (crit(7ri_i, K) f] Ur,,e f] B. 



Notice that lime_>oXe lies in lim£_5.o crit(7rj, Fe) and Ur^ C ext(C/, R). Similarly lim£_>oye lies in 
lime_j.o crit(7rj_i, V^e) and [/^ C ext(C/, R). Combining these with the fact that x^ and y^ are not 
infinitesimally close w.r.t e, we conclude that 

Qim crit(7ri, F^)) f^Ur, ^ (lim crit(7ri_i, K)) p| Ur, . 

The only if statement. Assume that 

(lim crit(7r„ F^)) f]U ^ (lim crit(7r,_i, F,)) f| U . 

We have to prove that there exists (x^, y^ Zs) E A^ and the lime_>o ||xe — y^H = lime_>o z^ > 0. 
First, we prove the existence of x^. Take x G ((lime_j.o crit(7rj, Ve)) fl ^ \ (li™e-5>o crit(7rj_i, Vs)) H U). 
By noticing that x G C/^, we deduce that 

X G ((lim crit(7r„ V^)) fl C/^ \ (lim crit(7ri_i, Ve)) fl uA . 

We claim that this implies that there exists Xg G crit(7rj, Ve) H Uri^e such that 



X = lim Xj G ( lim crit(7rj, Ve) ] (1 U, 



Indeed, since x G lim£_^o crit(7rj, y^) and / is non-negative over R(C)"i we deduce that there exists 
X£ G crit(7rj, Ve) n R(e)" with lime_>oXe = x (Lemma 18). 

Moreover, x G Vr^. But x G R(C) and the coefficients of the ^j's lie in R(C)5 thus 51 (x) > 
r/, . . . ,gs(x) > T]. If there exists and a /c, 1 < fc < s, such that gki^e) ^ V^ then this would imply 
that gkilinie^oXe) = 5'fc(x) < r]. This contradicts the fact that x £ Urj and so gk{^) > "H for all k. 
We conclude that x^ G C/^,£ and so x^ G crit(7rj, Ve) f] Ur^^e as we claimed. 

We also deduce that Xg ^ crit(7ri_i, K) H Ur^,e- Since in this case 



X = lim Xg G ( lim crit(7rj_i, Ve) ] (1 U 



which is a contradiction. Thus, there exists x^ G (crit(7rj, Fg) — crit(7rj_i, T^)) P| [/^^e P|i?e. The 
inclusion x^ G -Bg is a direct consequence of the fact that Be strictly contains the semi- algebraic 
set under consideration. 

Now we prove the existence of y^. Recall that x is obtained as the limit of x^ that it is 
not infinitesimally close, w.r.t e, to any point of the set crit(7rj_i, 14). By the assumption that 
crit(7rj_i, 14) Pi f/jj^g is not empty, there exists y^ G crit(7rj_i, y^) f] C/^^e that is not infinitesimally 
close to Xg, w.r.t. e. 

As there exist x^ and y^ that are not infinitesimally close w.r.t. e, there exists Ze G R(e) such 
that Ze > and ||xe — yep > z'^ and lime_s.o Ze > 0. 
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Complexity analysis. When DisjointPolar is called with i = 1 its run time is the one of 
LimitsOfOfCriticalPoints. By Lemma 16, Step 1 costs {s D)^""' arithmetic operations. 

Now we assume that i > 2. Steps 2-3 and Step 5 are symbolic manipulations which produce 
Boolean combinations of 0(s) polynomials of total degree 2D involving 0{n) variables. 

By Lemma 11, Step 4 costs at most [s D)^"'' arithmetic operations. 

The complexity of one-block quantifier elimination at Step 6 is (sZ?)'^^"' (Theorem 14). 

The output of the quantifier elimination procedure consists of (s-D) '"^ conjunctions of poly- 
nomials in Z[(^, 77, Z][e] of degree at most D^^'^> . Each conjunction involves at most {s D)^^'^> 
polynomials. 

After we apply UnivariateLimit (Step 8) we check if there is a realizable sign condition. According 
to Lemma 15 this costs (s D)^^' arithmetic operations. 

We call IsRealizable at most {s D)^^^' times (which is the number of conjunctions and thus the 
cardinality of \I'). Each call involves {s D)^"'' polynomials in Z[(^,r],Z] of degree at most Z)*^'") 
(Lemma 17). The arithmetic cost is (^s D^^""')^^^' = (s D)'-^^"'' (Sec. 5.2) which is also the cost for 
the whole for-loop and the algorithm. 

Statements on bit complexity when Z = Z are straightforward applying mutatis mutandis the 
same reasoning as above and using bit complexity results given in Theorem 14, and Lemmata 15, 
16, and 17. D 

5.3 Proofs of subroutines 

Proof of Lemma 11 . We decide the emptiness oiV {^U using the algorithm from [50, Proposi- 
tion 2.2]. The cost is (s S)^""' operations in Z. The Boolean complexity bound follows by combining 
[50, Proposition 2.2] with [13, Algorithm 13.1]. D 

Proof of Lemma 15 . The quantifier free formula of the input represents the following semi- 
algebraic set 

Se = {z^R{0{r^){e)\^} , 

where ^ is a Boolean formula whose atoms are polynomials in Z with coefficients in R{C,) {rf) {e) . 
We can express the limit of 5^ as e — )• 0, 5, using the language of the first order theory over the 
reals [13, Sec. 3.1], that is 

5 = lim5e = {z' G R(C)(?/)|(Vr >0)(3eo>0)(Ve)(Vz)zGcSA(0<e<eo^ \\z-zf < r^)} . 

In this way, to compute the limit it suffices to eliminate the quantifiers from the previous formula. 
For the elimination process, we treat ^ as a variable. To see that this is valid let <1> be the 
Boolean formula the describes z ^ S f\{Q < e < eq ^ \\z — z'W"^ < r"^). Then 

^^,^ = {(z,z',r,eo,e)GR(C>('?>'l^} • 

If we let C = t and r/ = /i in all the polynomials in $, then we get the set 

yl = {{(z,z',r,£o,£,t,/i)eR6|$}} . 

In this way 

7r_t(ext(A, R(C)(r?)) fji* = C,h = r^})= A^,, C R(C)(r/)^ , 

where 'K-th is the projection (z, z' , r, eq, e, t, h) 1— )• (z, z\ r, eq, e). 
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We notice that S is defined using a constant number of free variables, (<^, r], and Z'), a constant 
number of quantified variables, (z, r, eo, e), and a constant number of quantifier alternations. 
Following [13, Theorem 14.16], after we perform quantifier elimination, the output is a quantifier- 
free formula Vf=i A/Li hij t>ij 0, where hij e Z[C,r]][Z'] have degree at most S'^^'^\ I < (^(5)*^(^\ 
Ji < (£6)^^', and the complexity of the quantifier elimination procedure is {£6)^^' arithmetic 
operations in Z. 

If Z = Z and the input polynomials have coefficients of maximum bit size r, then the Boolean 
complexity is r {i5)^^^K D 

Proof of Lemma 16 . By assumption, crit(7ri,V^) is finite and H is non-negative over R(C)"- 
Thus, to compute (lime_5.o crit(7ri, Ve)) H f^ it is sufficient to 

(1) compute sample points in each connected component of crit(7ri, Ve) n R(C)(e)"; they will be 
encoded with a rational parametrization with coefficients in Z[^,e] 

q{T) = 0,Xi = qi{T)/qo{T), . . . , X„ = qn{T)/qo{T); 

(2) use this parametrization to compute their limits and make the intersection with U. 

The above parametrization is obtained using [13, Algorithm 13.1] with input H — e = ^^ = 
••• = ^^ = 0. Following mutatis mutandis the same reasoning as in [50, Proposition 2.2], we 
deduce that the arithmetic cost in Z of this step is (s(5)'^'"^ 

Step (2) is performed using with [13, Algorithm 11.20] (Removal of an infinitesimal) for compute 
the limit of the sample points as e — >■ and sign determination algorithms for univariate polynomials 
to obtain the intersection with U (see [13, Algorithm 10.13]). The overall complexity is (s(5) '"•' 
operations in Z [13, Chapters 10 and 11]. 

The Boolean complexity bound is straightforward from the above reasoning. D 

Proof of Lemma 11 . The subroutine is based on Algorithm 10.13 from [13]. If the degree of 
the polynomials is at most b and if there are at most i polynomials, then the total cost is ibr^^' 
arithmetic operations in Z[(^,r/] [13]. 

The operations that we need are computations of subresultant sequences for univariate poly- 
nomials in Z with coefficients in Z[^,7y]. The coefficients of the polynomials in the sequence have 
degree at most bP^^' in Q and r] [13, Proposition 8.49]. Therefore the complexity of the algorithm 
is also ^J*^'^-*, when we count operations in Z. 

The Boolean complexity is r ib^^> and it is due to the bit size of the integers that appear in 
the subresultant sequence. D 
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